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We describe a new method, based on an extension of the unitarity cutting rules proposed in Ref. [ ||], which 
is very efficient for the algorithmic evaluation of phase-space integrals for various differential distributions. As 
a first application, we compute the next-to-leading order normalized rapidity distribution of the CP-even and 
the CP-odd Higgs boson produced in hadron collisions through gluon fusion. We work in the heavy top-quark 
approximation; we find that the NLO corrections at the LHC are approximately 5% in the zero rapidity region. 


1. Introduction 

One of the major tasks of the Tevatron and the 
LHC is to discover and explore the so far inacces¬ 
sible Higgs boson sector of the Standard Model 
(SM). The discovery of a single CP conserving 
Higgs boson, as predicted by its minimal version, 
or a more prolific spectrum of Higgs bosons, char¬ 
acteristic to extensions of the SM such as the 
minimal supersymmetric SM (MSSM) or the two- 
Higgs-doublet-model (2HDM) will elucidate the 
nature of electroweak symmetry breaking. 

The dominant mechanism for the production 
of light Higgs bosons at hadron colliders is gluon 
fusion through a heavy quark loop. The pro¬ 
duction cross sections of both the CP-even (H) 
and the CP-odd ( A ) Higgs bosons are known ex¬ 
actly through next-to-leading order (NLO) in per¬ 
turbative QCD [ H and through next-to-next-to- 
leading order (NNLO) only in the infinite top- 
quark mass approximation [ 0 0. Note that, in 
the case of the CP-odd Higgs boson, this approx¬ 
imation is reliable for small values of tan (3 < 10, 
where the contributions of bottom-quark loops 
can be ignored. The double differential rapid¬ 
ity and px distribution for the SM Higgs boson 
has been calculated through NLO by means of 
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a fully differential Monte-Carlo program) 0, and 
analytically [ || in the case of non-zero pr- For 
the CP-odd Higgs the double differential rapidity 
and pt distribution was also calculated recently [ 

I- 

In this paper, we compute the NLO rapidity 
distributions for the production of the CP-even 
and CP-odd Higgs bosons analytically, including 
the virtual corrections at zero rapidity. This is 
formally one order lower in a s than the NLO con¬ 
tributions of Refs. [ 010 - (Some numerical re¬ 
sults for this distribution, including detector cuts, 
have been reported previously [ 0 ].) For the com¬ 
putation of the inclusive phase-space integrals for 
fixed Higgs boson rapidity we extend the method 
of Ref. [ g to accommodate the calculation of 
differential distributions. The idea is to replace 
the 5-function constraint on the phase-space by 
an “effective” propagator. This propagator de¬ 
pends on the constraint and in general differs 
from conventional particle propagators; however, 
if the constraint is polynomial in external mo¬ 
menta, the resulting Feynman integrals can ef¬ 
ficiently be dealt with by algebraic means. We 
will illustrate how this method works in the next 
Section. As a cross-check, we have also computed 
the rapidity distributions by explicitly integrating 
the finite remainders of the phase-space integrals 
after dipole subtraction [ 0. We found complete 



2 


agreement between the two methods. 


2. Method 


For the calculation we use the large top-quark 
mass (nit oo) approximation which is known 
to work extremely well [ even for relatively 
large Higgs boson masses. In this limit the inter¬ 
action of the Higgs boson with gluons is given by 
an effective Lagrangian [ 0-0 which is known to 
NNLO in the strong coupling constant. Keeping 
only the terms relevant to an NLO calculation, 
the CP-even and CP-odd effective Lagrangians 
read: 
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where G a is the gluon strength tensor, = 
€^ 1 ;a/ 3 G a P, H,A are the Higgs boson fields, and 
v ~ 246 GeV is the Higgs boson vacuum expecta- 
tion value. The Wilson coefficients C 1 ’ , defined 


in the MS scheme, are [ 0: 
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where a s = a s (p) is the MS strong coupling con¬ 
stant defined in the theory with n/ = 5 active 
flavors. 

We consider the collision of two hadrons with 
momenta Pi = Vs/2 (1, Ox, 1)) and P 2 = 
Vi/2 (1, 0 X ,-1), producing a Higgs boson h = 
{H, A} with momentum Ph = (P, prp,p 2 ). The 
rapidity Y is defined by: 
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The hadronic rapidity distribution is obtained 
from the partonic rapidity distributions by con- 
voluting them with appropriate parton densities: 
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The partonic rapidity distributions for the hard 
scattering of partons i, j with momenta pi = X\P\ 


and p 2 = x- 2 P 2 respectively, are obtained by inte¬ 
grating the hard scattering matrix elements over 
the phase-space of the final-state particles with 
the rapidity of the Higgs boson kept fixed: 
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In the center of mass frame of the colliding 
hadrons, the rapidity constraint can be written 
as: 

/ 2V . _ fi+fc \ = 2 y / (8) 

V E + pJ V Pi. Pi J 

with u = (xi/x 2 ) e~ 2y . 

At leading order in a s a sole Higgs boson is pro¬ 
duced; in this case momentum conservation ren¬ 
ders the phase-space integrals trivial. At NLO, 
the production of the Higgs boson is accompanied 
by a production of either a quark or a gluon. This 
makes the phase-space integrations more compli¬ 
cated, but they are still sufficiently simple to be 
done directly. (See for example Ref. [ |o| for the 
analogous case of Drell-Yan production.) How¬ 
ever, the brute force approach becomes cumber¬ 
some at NNLO and beyond. 

In this paper we show how to compute the 
NLO contributions using a method suitable for 
an algorithmic evaluation of the rapidity distri¬ 
bution at NNLO as well. The idea is to replace 
the (5-function constraint on the phase space in 
Eq. (|j) in terms of the imaginary part of an ef¬ 
fective “propagator”: 
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We can then map the constrained phase space 
integrals onto loop integrals in a manner sim¬ 
ilar to what was suggested for unconstrained 
phase-space integrals in Ref. [ [2). It is impor¬ 
tant that the constraint in Eq. (gj) is a polyno¬ 
mial in momenta; this property allows the ap¬ 
plication of multi-loop algebraic techniques, such 
as integration-by-parts and recurrence relations [ 
a to the integrals produced after the map¬ 
ping (|). 

At NLO, using Eqs. (|],^|), we can express all 
the phase-space integrals through linear combi- 













3 


nations of the following loop integrals: 


1 , . . . 

r d d k 1 

7 (2 ir) d A? ■ • ■ A? ' 

(10) 

where 



Ai = 

k 2 — m 2 ± id, 

(11) 
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(k + Pi +P 2) 2 ±iS, 

(13) 
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k ■ pi — uk ■ p 2 ± id. 
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The propagators Ai, A 3 and A 5 are “cut” 

accord- 


ing to Eq. (§). 


We now proceed to the reduction of the inte¬ 
grals of the above topology. It turns out that we 
can derive a sufficient set of recurrence relations 
through partial fractioning. Integration by parts 
is not needed for this calculation but it will be an 
essential tool at NNLO. 

The five propagators of the topology are lin¬ 
early dependent: 

Ai + A 3 A 2 — A 4 — s — T 7 ?, 2 , 

2 A 5 + A 3 ( 1 — u) + uA 4 — A 2 = (u— l)m 2 , 

where s = (pi + P 2) 2 ■ Using the above rela¬ 
tions we can eliminate both propagators A? and 
A 4 . It should be noted that partial fraction¬ 
ing produces terms with one or more of the cut 
propagators Ai,A 3 ,A 5 eliminated too. These 
terms have a zero contribution to the rapidity 
constrained phase-space integrals and we discard 
them. Finally, /(iq,..., 1 / 5 ) reduces to a single 
master integral 7(1, 0 , 1 , 0 , 1 ) = A'i in an algebraic 
fashion. Upon reinstating the <5 functions from 
the cut propagators, the master integral becomes 
the two-particle phase-space integral evaluated at 
fixed rapidity: 


X\ -> j [dh][d<j] ( 2 n) d S^ d \pi + p 2 - h - q) 
x 6 {h ■ [pi - up 2 ]) 

= [j/(l ~y)(l - z) 2 ] £ l-y + yz 
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The real radiation graphs are singular at z = 1 
and y = 0 or y = 1 . We extract the poles in e 
using identities of the form: 


,-i+e 
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Upon combining the real and virtual contribu¬ 
tions and performing the UV renormalization and 
mass factorization in the MS scheme, all the poles 
in e cancel and a finite result is obtained for the 
rapidity distribution. 


3. Partonic distributions 


We now present the analytic expressions for the 
partonic rapidity distributions. We write: 
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At leading order only the gluon-gluon production 
channel contributes: 

^g' AU0) = \ S(l-z)S(y[l- y ]). ( 22 ) 

At NLO we obtain contributions from the quark- 
antiquark, quark-gluon and gluon-gluon channels: 
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The above expressions are valid when the renor¬ 
malization and factorization scales are set equal 
to the mass of the produced Higgs boson. The 
full dependence of the partonic cross sections on 
those scales can easily be restored by solving the 
renormalization group and DGLAP equations us¬ 
ing the above expressions as boundary conditions. 

As expected, by integrating the partonic distri¬ 
butions over the rapidity Y, we obtain the par¬ 
tonic total cross sections calculated earlier [ [I|]. 


4. Numerical Results 

In this section we present numerical results for 
the NLO rapidity distributions of the CP-even 
and CP-odd Higgs bosons at the LHC and the 
Tevatron. We calculate the hadronic rapidity dis¬ 
tributions by convolving the partonic cross sec¬ 
tions of the previous section with the NLO parton 
distribution functions, as in Eq. (^[). The result¬ 
ing rapidity distributions, normalized to the total 
cross section, are shown in Fig. [I]. Fig. || shows 
the analogous plots for the Tevatron. It is clear 
from the plots that the corrections to the shape 
of the distributions are fairly small. For exam¬ 
ple, at zero rapidity, where the corrections are 
largest, they increase the LHC result by only 5%. 
For the Tevatron, the NLO rapidity distribution 
falls within the band of the LO distribution. The 
dependence on the factorization and renormaliza¬ 
tion scales is also small, suggesting that higher 
order perturbative corrections are unlikely to be 
large. 


This very stable behavior should be contrasted 
with the known fact that the NLO corrections 
to the total Higgs boson hadroproduction cross 
section are very large; they increase the cross sec¬ 
tion by approximately a factor 1.7. Our result in¬ 
dicates that since the shape of the distribution is 
very stable against higher order QCD corrections, 
it can be reliably predicted even by LO Monte 
Carlo event generators normalized to the NNLO 
results for the total cross section. This procedure 
should give a fairly accurate description of the 
Higgs rapidity distribution at the LHC. 

In the case of the CP-odd Higgs boson, the 
situation is rather similar. The partonic cross 
sections for the CP-even and the CP-odd Higgs 
bosons differ only in a single term proportional 
to 5(1 — z) with a small coefficient, Eq. (|2^); 
therefore the rapidity distributions for the CP- 
odd Higgs boson are numerically very similar to 
the distributions shown in Fig. 0 and FigJ|. 


pp -* H+X Rapidity distribution at LHC 



Yh 

Figure 1. CP even Higgs boson rapidity distri¬ 
bution at the LHC at leading (red) and next-to- 
leading order (blue) in perturbative QCD. 


5. Summary 

In this paper we computed the NLO rapidity 
distribution of CP-even and CP-odd Higgs bosons 
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pp -* H+X Rapidity distribution at Tevatron 



Figure 2. CP even Higgs boson rapidity distribu¬ 
tion at the Tevatron at leading (red) and next-to- 
leading order (blue) in perturbative QCD. 


produced at hadron colliders. We found that the 
NLO corrections change the rapidity distribution, 
normalized to the total hadronic cross section, 
only by a small amount. For example, at zero 
rapidity the NLO normalized distribution for the 
LHC increases by approximately 5% as compared 
to LO. The scale variation decreases by a factor 
of two, from LO to NLO. 

The phase-space integrations of the real radia¬ 
tion graphs with fixed rapidity of the Higgs boson 
are straightforward at this order in perturbation 
theory. However, traditional methods are very 
cumbersome for the evaluation of the Higgs bo¬ 
son rapidity distributions at NNLO. In this paper, 
we performed the first test of the extension of 
the method suggested in Ref. [ ||] for evaluating 
phase-space integrals using multiloop techniques, 
by applying it to the rapidity distribution of a 
hadroproduced color-singlet state. 

We are confident that the same method is 
tractable for the evaluation of the differential dis¬ 
tributions in more complicated cases, such as the 
rapidity distribution for Drell-Yan and Higgs bo¬ 
son hadroproduction at NNLO. This will be the 
subject of a future work. 
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